AiwtomlMr 


0 


Time-Dependent Simulations of Incompressible Flow in a 
Turbopump Using Ovreset Grid Approach 


Cetin Kiris 
Dochan Kwak 

NASA Ames Research Center 


European Congress on Computational Methods in Applied Sciences and 
Engineering (ECCOMAS 2001) 

Swansea, United Kingdom, September 4-7, 2001 


Outline 


tMsimntCulir 


► INTRODUCTION 

- Major Drivers of the Current Work 

- Objective 

► SOLUTION METHODS 

- Summary of Solver Development 

- Formulation / Approach 

- Current Challenges 

►PARALLEL IMPLEMENTATION 

► UNSTEADY TURBOPUMP FLOW 

► DISCUSSION 


Pag© 1 




A 

J'tomBssuttoCutor. 


Major Drivers of Current Work 


• To provide computational tools as an economical option for developing 
future space transportation systems (i.e. RLV subsystems development) 


Impact on component design 
Increase durability /safety 


Rapid turn-around of high-fidelity analysis 

Accurate quantification of flow 

(i.e. prediction of low-induced vibration) 


Impact on system performance => More complete systems analysis 

using high-fidelity tools 

Target 

Turbo-pump component analysis =* Entire sub-systems simulation 
Computing requirement is large: 

=>The goal is to achieve 1000 times speed up over what was possible in 1992 


Objectives 
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• To enhance incompressible flow simulation capability for developing 
aerospace vehicle components, especially, unsteady flow phenomena 
associated with high speed Turbo pump. 
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IN530 - Incompressible N-S Solver 
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** Parallel version : Based on INS3D-UP 

•MPI and MLP parallel versions 

• Structurel, overset grid orientation 

• Moving grid capability 

• Based on method of artificial compressibility 

• Both steady-state and time-accurate formulations 

• 3 rd and 5 th -order flux difference splitting for convective terms 

• Central differencing for viscous terms 

• One- and two-equations turbulence models 

» Several linear solvers : GMRES, GS line-relaxation, LU-SGS, 

GS point relaxation, ILU(0) 

•HISTORY 

** 1982-1987 Original version of INS3D - Kwak, Chang 
** 1988-1999 Three different versions were devoped : 
INS3D-UP / Rogers, Kiris, Kwak 
INS3D-LU / Yoon, Kwak 
INS3D-FS / Rosenfeld, Kiris, Kwak 


Time Accurate Formulation 
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Time-integration scheme 

Artificial Compressibility Formulation 

• Introduce a pseudo-time level and artificial compressibility 
* Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisf ied. 


Pressure Projection Method 

* Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 
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Current Challenges 
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• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important such as 

- Scalability (algorithmic <& implementation) 

- Portability, transparent coding etc. 

• Computing resources 

- ‘Grid* * computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using "super 
node" which is largely based on parallel architecture 


INS3D Parallelization 
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INS3D-MPI 

(coarse grain) 

T. Faulkner A J. Dccles 



• INS3D-MPI / Open MP *■ 

MPI (coarse grain) ♦ OpenMP (fine grain) o 

Implemented using CAPO/CAPT tools * 

H. Jin AC. Kins 


OpenMPl 
L, threads I 


Group 1 
MLPftocwsf 


Group N 
ULP Process 2 


INS3D-MLP 
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Time (sec) per iteration 
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SSME-rigl /Initial start 
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TIME STEP 96 / Impeller rotated 42-degrees 


VELOCITY MAGNITUDE 


- 34.3 Million Points 

- 800 physical time steps in one 
rotation. 

*One physical time-step requires 
less then 20 minutes wall time 
with 80 CPU's on 0rigin2000. 
One complete rotation requires 
one-week wall time with 80 CPUs. 
^Currently I/O is through one 
processor. Timing will be improved 
with parallel I/O since time- 
accurate computations are I/O 
intensive. With futher 
improvements several impeller 
rotations can be completed in one 
week 


Summary 
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•Unsteady SSME-rigl start-up procedure from the pump at rest has been 
initiated by using 34.3 Million grid points. 

• Computational model for the SSME-rigl is completed. Moving boundary 
capability is obtained by using DCF module in OVERFLOW-D. 

• MPI /Open MP hybrid parallel code has been benchmarked. 

• MLP shared memory parallelism has been implemented in INS3D, and 
benchmarked. 

• MLP/OpenMP version requires 19-25% less computer time than 
MPI/OpenMP version. Pin-to-node for MLP version is implemented. 40% less 
computer time is required in the new version. 

• Time-accurate features of methods designed for 3-D applications are 
evaluated. An efficient solution procedure is obtained. 

• Work currently underway 

•Unsteady SSME-rigl simulations by using 34.3 Million grid points. 

• Experimental measurements at NASA-MSFC. 
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